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Abstract: Battery is a major hardware component of wireless sensor networks. Most of them 
have no power supply and are generally deployed for a long time. Researches have been done on 
battery physical model and their adaptation for sensors. We present an implementation on a real 
sensor operating system and how architectural constraints have been assumed. Experiments have 
been made in order to test the impact of some parameter, as the application throughput, on the 
battery lifetime. 

Key-words: battery, wireless sensors network, Contiki, Cooja 



* Universite de Lorraine 



RESEARCH CENTRE 
NANCY - GRAND EST 

615 rue du Jardin Botanique 
CS20101 

54603 Villers-les-Nancy Cedex 



Estimation en ligne de la Duree de Vie des Batteries pour 
les Reseaux de Capteurs sans Fil 

Resume : La batterie est un composant materiel majeur dans les reseaux de capteurs sans 
fil. Des recherches ont ete menees sur les modeles physique des batteries et leur adaptation 
pour les capteurs. Nous presentons une implantation sur un veritable systeme d'exploitation de 
capteurs et comment les contraintes d'architecture ont pu etre pris en compte. Des experiences 
ont ensuite ete menees pour tester l'influence de certains parametres comme le debit d'application 
sur la duree de vie des batteries 

Mots-cles : batterie. reseaux de capteurs sans fil, Contiki, Cooja 
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1 Introduction 

Wireless Sensors Networks have a growing place in industrial and domestical domains. Gener- 
ally sensors need to be small, cheap and easy to deploy on any physical environment. These 
constrained networks are made possible by the processor size reduction and high performance 
battery. Nevertheless these networks have to be care of their energy consumption because they 
are usually deployed for a long time or because sensors are so small that their battery has a 
limited capacity. The growing interest on battery lifetime estimation for WSN appears in the 
standardization work on these network [B]. They defines a routing protocol with one metric being 
the remaining battery level |20j . 

Battery models that relate its use with its lifetime have been proposed since embedded systems 
(laptop, cellular phone) are widespread. All of them take into account two phenomenons occuring 
in battery cell; the Rate Capacity Effect and the Recovery Effect [S]. The first gives the energy 
consumed under a constant current load, when transmitting for example, and the second gives 
the energy recovered during inactivity or low current load. At a constant current load, oxidation 
at the anode induces reduction at the cathode. The reduction decreases the concentration of 
positive ions near the cathode and so the available energy. But during inactivity or low current 
load the positive ions near the anode have time to move toward the cathode, thus increasing 
the available energy and so the battery lifetime. A stochastic model is given in [5], but the 
more accurate modeling is given by the analytical equation ((T|) [T31 HH Qj)] . a denotes the total 
battery capacity that is equals to the current load i(r) consumed since the use of the battery 
until the time L (the battery Lifetime) and the current that was unavailable because positive 
ion concentration was not sufficient at the cathode. The (3 parameter is the diffusion coefficient 
of electroactive species inside the battery. 

!(r)dr + 2V / i{T)e~^ m ^ L - T) dT (1) 

Use of this model can be found in Q2J [TS] to schedule tasks for embedded system or to 
simulation [19J . It is widely used for example by a circuit-based battery model or in WSN 
routing simulation [17] . The model can not be implemented as is on WSN nodes albeit some 
results achieved using the first part of equation ([T]) [3_J [7] . In so doing, they do not take into 
account the recovery process that occurs during idle time despite it is generally up to 90% of the 
WSN node timelife. 

The contribution of [12 is an approximation of the equation (TTJ) that can be implemented on 
WSN node. Its global equation © recursively computes a(L n ) that is the consumed charge in 
mA-min at time L n . This computation depends on the consumed charge at time £ n _i, the time 
interval A = L n — L„_i being constant. 

n / n—1 \ 

a(L n ) = ^2 I k 6k + A I <r(L n _i) - ^ I k 5 k + 2I n A(L n ,L n _ 1 + <J fc ,L n _i) (2) 

k=l V fe=l / 

With this approach, one can know the remaining charge in the battery at time L n by the difference 
between a(L n ) and the initial battery parameter a of the equation (TTJ). The recovery effect of 
the battery is computed through the function A of the equation @ that is approximated by the 
use of an / function: 

A[L n , L„_i + d k> L n _!) - 2^ Q2^2 - ~Q2 ff2~ W 

rn—1 
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This approximation depends on the idle time of the mote during the A interval, given by the 
difference v = A — 5k (y is the idle time of the L„_i period and 5k its activity time). Following 
|12) , the / function can be also approximated by the use of the equations (QJ to © (see [TDl [TT] 
for details): 



(5) 



2^ 



/(A) _ ~ 

771=1 

The /3 parameter (physical diffusion coefficient) in the equations Q and (J5]) is a constant 
value that depends on the battery and is computed from several tests of discharge and the least 
squares estimation method (along with a). The approximation of the equation ((5|) concerns y/v 
by a mean value : a, computed offline. In our case the mean value a is fixed at 90% of A, that is 
the value we measure on several tests. The last approximation of the equation ^ is computed 
with the m parameter ranging from 1 to 10 because the sum quickly converges. 

The recursivity of the model depends on the A parameter defined as the ratio ~^r£ ± ^~gp f° r 

each n but that can be bounded by the value e~P A computed off line. 

Our contribution in this paper is to fill the last gap between the approximated model and 
its implementation within an existing operating system for sensors. We use physical nodes to 
validate our energy consumption retrieval process and observe in simulated world how battery 
reacts face to network events and configuration parameters. The remainder of this paper is 
organized as follow. The next section is focalized on the elementary operations of WSN node 
and their current draw. Details of the implementation are given in the section [3] and first results 
are shown in the last part. We give some conclusions and plan future works at the end of the 
document. 



2 Current draw 
2.1 Linear draw 

The first term in equation ([2J is the sum of products between current load Ik and time interval 
5k since the begining of the battery livetime (when k = 1) until the time L n . The 5k time interval 
is a sub interval of A during which the battery was used by the node (5k < A). To compute 
this value we have to know which components of the motes use the current along the time and 
at which current rate. 

Mote's current information can be obtained from datasheet document. For example, the table 
Q] shows the current load for two mote types (Sky or WSN430 motherboard, respectively with 
CC2420 and CC1100 communication chips). The columns are the usual states of the duty cycle 
for any mote : CPU, LPM, TX and RX respectively for processor, low power mode, transmitting 
and receiving. Currents are given in milli- Ampere (mi). The given values for TX and RX are 
linked to signal power and throughput configuration of the mote. 

The left part of figure [1] shows a short time interval of a real WSN430 mote with the Contiki 
operating system [5] on the senslab testbed [2]. We plot these states by logging start and stop 
time on mote serial output that is redirected by senslab to a TCP connexion. These times are 
accessibles by an energy related application programming interface on Contiki. The CPU and 
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Mote 


CPU 


LPM 


TX 


RX 


Sky 


1.8 


0.0545 


17.4 


18.8 


Wsn430 


2 


0.02 


16.1 


15.2 



Table 1: Motes current draw 



LPM are not overlapping themselves and fill all the time duration. These states are related to 
the mainboard activities and so either the CPU is active or the low power mode is raised. TX 
and RX states are related to the radio activities and they are neither used every time nor they 
overlap each otther. One can note that TX mode only occurs during CPU and that RX mode 
begins during CPU but can spread on LPM. 
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Figure 1: Mote modes and currents 

The senslab testbed also provides an external polling of the current consumed by a mote 
during an experiment. The right part of figure [T] shows the curves of measured current at the 
mote endpoints and the duties cycles as in the left part (we just leave TX and RX mode that 
have the main current values). One can see that start and stop times of duties match well with 
measured current. 

Consequently, the current draw of a A = 5c pu + &lpm intervall is the result of the equation 
10 where each C s t a te is the current load (cf.table [T]) and each S s t a te is the sum of all state 
sub-intervals inside A. 

IkSk = Ccpu-dcpu + Clpm-Slpm + Ctx-Stx + Cpx-Spx (7) 

3 Implementation constraints 

Even with the help of |12| . we must be care of mote limitations on code size and arithmetics 
capabilities. The strong constraint is that MSP430 can not handle floating point values but only 
signed or unsigned integer. 
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3.1 Offline computation 

Approximations given make use of y/n or ir 2 and we had to round them with a new unit as it is 
show on tabled Most of them are just a factor of thousand to remove floating point with saving 
enougth precision. The computed values v, cq and A are computed with a A of 2s. 
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9.869 


1.772 
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1.337 


0.967 


0.03 


0.173 


2.886 
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9869 


1772 


10 


1337 


967 


30 


173 


2886 



Table 2: Real to Integer values 



3.2 Online computation 

We compute the consumed current every A time interval. The values 8cpu,lpm,tx,rx are 
provided in milli second and we use the current load of the table [T] with a factor of thousand. 



We compute v with the low power mode period and the radio use : 

Slpm - (Stx + Srx) if Slpm > (Stx + Srx) 
else 



(8) 



If the radio is more used than the CPU is idle then there is no idle time for the battery. 

The varying part of the recovery function A (equation [5]) is rewritten in order to lose as little 
precision as possible : 

f(v) 10 4 /3 2 i/ + 7r 2 2.10 3 -i2/V^ 

/3 2 /3 2 12.10 5 [ ' 

At this step we had to use a 64 bits intermediary value to get the result (MSP430 compiler allows 
such data type but they are actually emulated with several 32 bits values). 

Finaly, the remaining energy is computed with the values above and the preceding remaining 
energy value (a A before). We compute the ratio with 255 as 100%. This value is given in 
RPL routing metric recommandation |20| and we use the present work to implement and test 
an energy-based routing plane in an other work beyond the scope of this paper. However we 
compute the remaining energy with a precision of 5 numbers, that is from 255. 10 5 , so one can 
observe energy consumption at a very fine coarse. 

The code size of the presented implementation is about 12Kb on the 50Kb allowed by our 
processor. At the memory level, we save on memory six 32 bits values from one computation to 
the other (times of CPU, LPM, TX, RX and previous values for <t(L„_i) and X)fc=i IkSk)- 



4 Simultation Results 

For our simulation we use the Cooja WSN simulator [5] and test a network of nodes build 
on the sky mote platform. These nodes use IPv6 and the RPL routing protocol to perform a 
simple collect application during which one or more nodes, the Senders, periodically send sensing 
informations to a receiving node called the Sink, configured to be the root of the RPL routing 
tree. All node have an initial battery charge of 880mAh. Senders send one data packet every 
second and use the contikiMac radio duty cycle We use a linux Ubuntu 11.0 desktop computer 
with a 3GHz Intel processor and AGb of memory. 
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4.1 Booting node 

We start our experiment with the observation of the node booting process. The figure [5] puts 
together the evolution of the battery lifetime and the time spends by activities of the node. 
During the first minute, on the part (a), one can note a large decrease of the batter}0to correlate 
with the strong activity on the part (b) for the same period. Mainly this activity is related to 
system and network initialization (remains the strong power for TX and RX of the table [T]). 
The recovery of battery is visible at this very close level, during the second minute and the two 




0123456 0123456 
Time (mn) Time (mn) 



(a) (b) 
Figure 2: First minutes of battery lifetime 

following ones. These recovery periods match with less activities period. 
4.2 Network building 

We continue the simulation following the figure [3] (a) where after ten minutes of simulated time, 
we add five nodes "under" the non-root node to observe the cost of the routing tree building 
and traffic cumulation. The node 1 is the Sink and all other are Sender. We leave these nodes 
sending and forwarding (only for the node 2) data packets during ten minutes and remove the 
five nodes previously added. This can occurs in real life for example with mobile nodes or with 
perturbed communication environment. A last point of interest, at the thirty fifth minute, is 
when a node lost its destination node (its default next hop router) because the network has to 
build itself again. 

The figure [31(b) plots the remaining battery lifetime of the node 2 along our experiment. The 
recovery curve is our implementation and the linear is a simple additive function of current 
draw. The first strong decrease of the tenth minute is mainly due to traffic overhead and very 
few is from RPL control messages. Once the routing plane is established, between the tenth and 
twentieth minute, the battery decreases with a greater slope and with more micro variations. 
When removing nodes at the twentieth minute, the node 2 shows a quick and strong recovery 
of its battery lifetime. As given by the model, the fall of current draw lets the battery recovers 
some of its current charge. During the fiftheen following minutes the node sends one packet per 
second to the Sink and the curve is similar to the begining (without booting process) but at a 
lower battery level. 

Hhat is only 0.0005% of the total 
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Figure 3: Battery lifetime and RPL 



At the thirty fifth minute we remove the parent of the node 2. The very large decrease of 
battery is again mainly related to a strong use of communication. The node has no more IPv6 
routeur neighbor then it polls the network with neighbour discovery messages and after three 
minutes, the node decides to build a new network and just send few RPL control message. One 
can there observe a very significant recovery at the part (b). 

The linear model is very less reactive and never grows, as expected. This strong difference 
suggest it could over or under estimate the lifetime depending on how the other sensors behave. 

4.3 Implementation Accuracy 

The approximation we made on basic values and the rounding effect of divide operations lead to 
residual errors all along the lifetime estimation. We have implemented a floating point version 
of the model with a perl script and compare its result with the computation made by the mote. 
The figure [4] shows a snapshot of a WSN like those of the figure [3]a with the two battery lifetime 
estimations. The mote computes values a bit lower than floating point but the difference stays 
constant. Strong variations (between the minutes 30 and 35) have a much more amplitude but 
do not deviate the lifetime. These results enforce our confidence on our implementation as it 
follows the theory of battery consumption and especially the recovery one. We continue research 
on the simulation for long time duration. 

4.4 Network life 

Estimation of sensor lifetime is a prior information for their deployement. Sensors networks can 
be in place for several years before be replaced or recharged and so battery should be chosen to 
match the expected duration. 

The Contiki operating system provides some Radio Duty Cycle (RDC) implementations we 
have used in order to compare their energy consumption. The part (a) of the figure [5] shows 
these experiments. All these RDC are asynchronous and packet oriented. 

• the contikiMAC [3] allows the greatest battery lifetime. It can sleep up to 99% of the 
time and was measured as been ten times less energy consumer than the following X-MAC. 
When a node has to send, it transmits the packet several times, until the receiver return 
an acknowledge. 
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Figure 4: Floating point vs Integer battery lifetime computation 



• X-MAC [TJ is a more consumer; the sender sends several preambule packets until acknowl- 
edge reception and then sends the data packet. 

• CX-MAC (Compatibility X-MAC) is a variation of X-MAC that is provided by Contiki to 
be usable on more radio chips. 

• sicslowmac (also in Contiki) handles 802.15.4 frames and has its radio alway open. Note 
that with a TX at ~ 20mA and a battery at 880m Ah, the theory gives ^ = 44/i (cf table 
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Figure 5: Network life 



We have also measured the throughput effect on battery lifetime and plot them in part (b) of 
the figure [5] The contikiMAC RDC was used because it is the less battery consumer and so the 
impact of the throughputs is more discernible. Numbers within the part (b) are the throughput 
of the application, from 60 packets by minute (pkts/mn) to 1. Each packet has a size of 87 
bytes and the application is parameterised by a number of seconds between two packets. The 
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figure shows that the lifetime is roughsly linear with the throughput but other tests with less 
throughput than Ipkt/mn have no significantly improved the lifetime. The network traffic is then 
negligeable against sensor base activity (indeed the routing protocol is the main radio user). 

These figures are very close to a line (only a zoom shows their are not) and we use a linear 
regression with the least square estimation technique to envision the time when the battery will 
be completely empty (i.e. with 0% of remaining energy). The table [3] contains the battery 
lifetime estimation for each lines in figure [SI The short lifetimes of these networks (no more than 
four months) is related to the light battery we have used. Battery of node like Sky are usually 
around 2000mAh and should allow WSN be operational for one year. Nevertheless real batteries 
have a secure cut off level that prevent them to be fully discharged because they will be damaged 
and not able to be charged again. 



RDC 


Battery lifetime 
(days) 


Throughput 

(pkts/mn) 


Battery lifetime 
(days) 


contikiMAC 


128 


60 


77 


X-MAC 


30 


20 


105 


CX-MAC 


26 


12 


113 


sicslowMAC 


1.8 


2 


124 


Throughput : lpkt/mn 


RDC : contikiMAC 



Table 3: Battery life estimation 



5 Conclusions 

This paper presents a first implementation of a battery lifetime estimation inside an existing 
operating system for WSN. Based on a recognized theoritical battery model, our work has show 
how sensor internal architecture impacts on concrete realization. Our simulations help us to 
better understand how sensor use their battery in bootstrap process or during networking. We 
believe this work is usefull for many applications as routing optimization that is a work we 
currently plan. We also have to launch long time tests establishing a relation between the 
reamaining energy and the cut off level of a battery. Moreover several battery technologies have 
to be tested to enforce our implementation model. 
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